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ABSTRACT 

We report the discovery of a substantial stellar overdensity in the periphery of the 
Large Magellanic Cloud (LMC), found using public imaging from the first year of the 
Dark Energy Survey. The structure appears to emanate from the edge of the outer 
LMC disk at a radius ~ 13.5° due north of its centre, and stretches more than 10 kpc 
towards the east. It is roughly 1.5 kpc wide and has an integrated E-band luminosity 
of at least My = —7.4. The stellar populations in the feature are indistinguishable 
from those in the outer LMC disk. We attempt to quantify the geometry of the outer 
disk using simple planar models, and find that only a disk with mild intrinsic ellipticity 
can simultaneously explain the observed stellar density on the sky and the azimuthal 
line-of-sight distance profile. We also see possible non-planar behaviour in the outer 
disk that may reflect a warp and/or flare, as well as deviations that resemble a ring-like 
structure between ~ 9° — 12° from the LMC centre. Based on all these observations, 
we conclude that our remote, stream-like feature is likely comprised of material that 
has been stripped from the outskirts of the LMC disk, although we cannot rule out 
that it represents a transient overdensity in the disk itself. We conduct a simple N- 
body simulation to show that either type of structure could plausibly arise due to the 
tidal force of the Milky Way; however we also recognize that a recent close interaction 
between the LMC and the SMC may be the source of the stripping or perturbation. 
Finally, we observe evidence for extremely diffuse LMC populations extending to radii 
of ^ 18.5 kpc in the disk plane (« 20° on the sky), corroborating previous spectroscopic 
detections at comparable distances. 
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1 INTRODUCTION 

The Large and Small Magellanic Clouds (LMC and SMC) 
are the nearest examples of star-forming, gas-rich dwarf 
galaxies. They sit well within the halo of the Milky Way, 
at Galactocentric distances of ~ 50 kpc and « 60 kpc, re¬ 
spectively. Their proximity renders them ideal targets for 
investigating a huge variety of different astrophysical prob¬ 
lems; their stellar populations can be easily resolved with 
ground-based facilities, and studied in exquisite detail. 

A critical ingredient in understanding many of the 
peculiar features of the Magellanic System - for ex¬ 
ample the striking off-centre bar in the LMC (e.g. 
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Ide Vaucouleurs &; Freemanl Il972l : Ivan der Marell l200lh , the 
broad age gap in the LMC star cluster system (e.g., 
|Da Costal 11991 L the vigorous ongoing star formation in 
both galaxies, the irregular morphology of the SMC (e.g., 
IWesterlundlIl997 ), and the Magellanic Stream and Bridge 
of Hi gas (e.g.. iMathewson et al ] |l974 IPutman et alJll99a . 
120031 : iNidever et al.ll2008l . 1201(f )- is the interaction history 
of the LMC and SMC, both with each other and with the 
Milky Way. This is still not well constrained. It has only 
been during the last decade, with precision proper motion 
measurments from the Hubble Space Telescope that evidence 
has mounted that the Clouds have only recently fallen deep 
into the gravitational pote ntial of the Milky Way (e.g., 
iKallivavalil et all 120061 . l2013h , and may indeed be on their 
first close passage. As a result it is possible that close en- 
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counters between the LMC and SMC have been more in¬ 
fluential in shap i ng their evolution and appearance (e.g., 
iBesla et al.ll2012l : iDiaz fc Bckkill2012h . 

The outskirts of the LMC and SMC are of particular 
interest. These low-density regions are most susceptible to 
disturbances due to interactions between the two galaxies 
and due to the tidal force of the Milky Way. Key questions 
are the extent to which the structure of the outer LMC 
disk exhibits evidence of distortion, the nature and extent 
of any stellar component stretching between the Clouds in 
the Magellanic Bridge region, whether there is a stellar com¬ 
ponent associated with the Magellanic Stream, and whether 
either of the Clouds possesses a spheroidal halo component 
or shows evidence for the presence of disrupted satellites. 
The latter may be expected because in the ACDM cosmolog¬ 
ical model structures are formed hierarchically at all scales, 
such that a number of the dwarfs in the Local Group could 
have experienced significant m erger events within the past 
several Gyr (e.g., iDeason et al.ll2014h . 

Studying the remote peripheries of the LMC and SMC 
is observationally extremely challenging, requiring excep¬ 
tionally low surface brightness features (V > 30 mag 
arcsec -2 ) to be traced over large swathes of sky. Although 
early studies utilising photographic plates were able to ex¬ 
plore th e morphology of the Clouds a t intermediat e radii 
(e.g., Ide Vaucouleurs & Freemanl Il972l : llrwinl Il99lh . it is 
only quite recently, with the advent of wide-held imagers 
and multi-object spectrographs mounted on medium and 
large aperture telescopes in the southern hemisphere, that 
the problem of tracing the LMC and SMC to their outer 
limits has become tractable. 

The structure of the LMC is generally accepted to 
be that of a planar disk viewed at an angle inclined 
to the plane of the sky. Pioneering work on the shape, 
structure and extent of the LM C disk was conducted by 
Ivan der Marel &; Cionil (|2001 j and Ivan der Mare! (l2Q0lh us¬ 


ing newl y-available contiguous near-infrared data sets from 
2MA SS (jSkrutskie et al.l 2006h and DENIS (|Epchtein et al.l 
Il997h . that provided, for the first time, the opportunity to 
map the outer regions of the LMC via resolved red giant 
branch (RGB) and asymptotic giant branch (AGB) star 
counts - i.e., using (old) stellar populations characteristic 
of the LMC periphery. These studies showed that the outer 
disk, or at least its old component, appears rather smooth 
all the way to the edge of the region studied at ~ 7° from 
the LMC centre, and found a best-fit inclination angle of 
i = 34.7° =b 6.2° and exponential scale length of « 1.3 — 1.5 
kpc. These values are consistent with those obtained by 
other studies conducted before and since at more c entra l 
radii using similar t r acers (e.g., Weinberg k, Nikolaevl l200ll : 
Idsen fe Salvkll2002l : iRubele et aLl2012f T as well as a variety 
of alternatives (such as Cepheids and Hi gas), and incor¬ 
porating both photometri c and kine matic techni ques (see 
the exten sive discussio n i n Ivan der Marel &; Cionil l200l| j. In 
addition, Ivan der Marel fc Cionil (l200lh observed tentative 
evidence for a decrease in the inclination angle with galacto- 
centric radius, together with a twisting of the position angle 
of the line of nodes - both possibly in dicative of war ping of 
the LMC disk at large radii. Moreover. Ivan der Marel (1200il l 
found that the LMC disk likely has an intrinsic elongation, 
perhaps induced by t he tidal for c e of th e Milky Way. 

More recently IBesla et~ahl ( 20161 1 have presented a 


panoramic view of the LMC in diffuse (unresolved) light out 
to ~ 8° using long exposures from small telescopes. They 
found evidence for faint stellar arcs and multiple spiral arms 
in the northern part of the disk, but with no comparable 
counterparts in the south. Their numerical models suggest 
this type of asymmetric arrangement can arise as the result 
of repeated close interactions with the SMC. 

The extent and structure o£ th e LMC beyond the 
edge of the reg i on stu died bvlva n der Ma rel &; Cionil (|200 il l. 
Ivan der Marelf (| 200 il l , and IBesla et al.l (|2016l l at ~ 7 — 8° is 
not well explored or understood. On the one hand, spectro¬ 
scopic measurements covering a significant number of dis¬ 
joint fields at a wide variety of azimuthal angles and ra¬ 
dial distances report kinematic detection o f LMC RGB stars 
out to ~ 20° (jMaiewski et al.lIl999l . l2009l l: moreover, LMC 
populations have also been clearly detected, both kinemat¬ 
ically and photometrically, in the foreground of the Carina 
dwarf at a radial distance of ^ 22° from the LMC centre 
dMaiewski et al. 1l200c l: iMunoz et al.ll2006l : iMcMonigal et a l.' 
12014 l. lMaiewski et all (|2009h hypothesise that these remote 
stars are indicative of a classical pressure-supported stellar 
halo surrounding the LMC. On the other hand, deep pho¬ 
tometric studies tracing old main s equen ce turn-off stars 
in stripes to the nor th (|Saha et al.1 l2010h and north-east 
([Balbinot et al.l l2015h observe a clear exponential decrease 
in surface density out to ~ 15° — 16°, suggesting that the 
smooth LMC disk extends at least 10 scale-lengths from 
its centre. Beyond this radius, neither work observes defini¬ 
tive evidence for LMC populations, and each posits that the 
LMC must be truncated , perhap s t idally, at ap proximately 
this point. Most recently. iBelokurov fc Koposov (2016) have 
probed to even larger radii in the north and east using blue 
horizontal branch stars, and report the discovery of an ap¬ 
parently rather substructured distribution reminiscent of a 
stellar halo extending to at least ~ 30°. 

In this paper we present a detailed investigation of the 
northern outskirts of the LMC, utilising images taken during 
the first year of the Dark Energy Survey (DES) and released 
publically through the NO AO science archive. In particu¬ 
lar, we explore the geometry and extent of the LMC disk 
and periphery, and report on the discovery of a substantial 
stream-like stellar substructure situated to the north and 
north-east of the LMC, and stretching to a radius of at least 
20° from the galaxy’s centre. 


2 DATA 

Our analysis is based on photometry derived from publically 
availa ble i mages o bt ained in the first year of DES operations. 
DES ([Abbott et al.l l2005l j u tilises t he 52 0 megapixel Dark 
Energy Camera fPECarm lFlaueher et al.ll2015T ) mounted on 
the 4m Blanco Telescope at the Cerro Tololo Inter-American 
Observatory in Chile. The DECam imager consists of 62 
CCDs arranged into a roughly hexagonal mosaic spanning 
a total field of view of « 3 deg 2 with a pixel scale of 0.26" 
pixel -1 . The DES survey area will ultimately spread over 
~ 5000 deg 2 of the southern sky in grizY\ the first-year 
imaging spans around ~ 2200 deg 2 with mostly single-epoch 
coverage. This includes a region sitting predominantly to the 
north and west of the LMC, covering position angles between 
~ —90° and ~ 45° at radii typically outside 7.5° but reach- 
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ing as close as ~ 6° to the LMC centre. Overall, roughly 
three-eighths of the LMC disk beyond 7.5° is covered by the 
year one DES imaging. 

We utilise the photometric catalogue described by 
IKoposov et al.1 ( 2015 b and refer the interested reader to 
that paper for a full description of its construction. In 
brief, source detection and measurement across all images 
was c onducted using the SExtractor and P SFEx software 
(e.g., iBertin k, Arnoutsl Il996l : iBertinl l201lh . The individ¬ 
ual detection lists were calibrated to the SDSS scale using 
DR7 of the APASS survey, cross-matched to remove dupli¬ 
cations, and merged across different pass-bands into a fi¬ 
nal catalogue. Star-galaxy separation was undertaken using 
the SPREAD_model and SPREADERR_model para mete rs pro¬ 
vided by SExtractor (as defined in lDesai et al.ll2012l ). Gen¬ 
eral improve ments were made to the processing over that 
reported by IKoposov et a D (l2015h . including an enhanced 
photom etric calibration utili sing a n “uber-calibration” pro¬ 
cedure ([Padmanabhan et al.ll2008h , extension of the process¬ 
ing to incorporate z-band data, and the use of all DES imag¬ 
ing from the 2013-2014 season (DES Yl) resulting in some¬ 
what increased areal coverage. 

For the present analysis, only objects with clean stel¬ 
lar detections in the g , r, and i filters were retained. This 
results in slightly brighter com pleteness limits than those 
reported bv IKoposov et al.1 (201 5) - by examining the posi¬ 
tion of the peak in source number counts we estimate that 
the catalogue starts to become significantly affected by in¬ 
completeness around g ~ 23.0, r ~ 22.9, and i ~ 22.2. 

For all sou rces, for eground red dening was corre cted us¬ 
ing the maps of Schlegel, Fmkbeiner k Davis (1998) and the 
extinction coefficients from lSchlaflv k Finkbeinen ([201 lh ap¬ 
propriate for the SDSS photometric system - i.e., 3.303, 
2.285, and 1.698 for the g , r, and i bands respectively. 
In the direction of the northern outskirts of the LMC, 
the reddening is very mild with a median colour excess of 
E(B — V) « 0.05. 


3 ANALYSIS 

3.1 Color-magnitude diagrams 

In the leftmost panel of Figure [I] we show a Hess diagram for 
all stars in our catalogue lying in the range 9° — 13° from 
the centre of the LMC. Most prominent is the old stellar 
population of the outer disk, which is well described by an 
isochrone f rom the Dartmouth Stellar Evolution Database 
([Potter et al.l l2008h of age 11 Gyr, [Fe/H] = —1.1, and 
[a/Fe] = 0.0. Our photometry reaches ~ 1 magnitude be¬ 
low the main-sequence turn-off belonging to this popula¬ 
tion. A distinct red clump is visible (marked with a box), 
but only a very faint blue horizontal branch, indicating that 
the peak of the metallicity distribution function sits around 
[Fe/H] % -1 or higher. This is consistent with pr evious 
investigations of ou ter LMC fields (e.g., ISaha et al.l l2010l : 
I Carrera et al.l 1201 lh . A trace younger population may ex¬ 
tend to the blue but this is completely sub-dominant. The 
rightmost panel of Figure [I] shows a Hess diagram for all 
stars in the range 7° — 9° from the centre of the LMC, 
and a substantial younger component is clearly present. This 
demonstrates that such populations are almost completely 


confined to a radii smaller than ^9°, at least to the north 
of the L MC - again consistent with previous studies. For 
example, iBalbinot et a D d2015h observed that populations 
younger than ~ 4 Gyr were truncated at ~ 9° or 10° from 
the LMC centre in their survey fields to the north-west. 

Returning to the leftmost panel of Figure [l] the distance 
modulus required to shift the fiducial isochrone onto the 
ridgeline of the dominant old stellar population is ~ 18.35. 
Taken at face value, this suggests that the region of the LMC 
imaged in the year one DES footprint mostly lies somewhat 
closer than does the galaxy centre, which has a distance 
modulus of ~ 18.5 (throughout this work we ass ume the dis¬ 
tance to the LMC centre to be Do = 49.9 kpc. lde Griis et ahl 
l20l4) . It is important to note, however, that our measure¬ 
ment is subject to systematic uncertainty of at least 0.1 mag, 
due to the fact that we do not know precisely the age, metal¬ 
licity, and a-element abundance of the stellar population(s) 
in the CMD, and because distances derived using isochrones 
can be rather model-dependent. 

Another notable feature of our CMD is that both the 
sub-giant branch and the red clump of the dominant old 
population are vertically quite broad. While this could be 
due to the presence of stars with a mix of metallicities and 
(old) ages, it might also reflect a spread in the distances to 
LMC stars in our catalogue. Indeed, we find the latter possi¬ 
bility to be the most likely explanation - in the two central 
panels of Figure [I] we plot Hess diagrams for stars lying in 
the radial range 9° — 13°, but at position angles ^ 25° (cen¬ 
tre) and ^ —55° (right) - i.e., to the extreme east and west 
of the imaged region. A significant distance gradient is ap¬ 
parent, in the sense that stars further to the east are closer 
than those to the west. The sub-giant branch in the middle 
panel lies ~ 0.15 magnitudes above the fidicual isochrone, 
while that in the right-hand panel sits ~ 0.1 mag below the 
fiducial line; there is also a striking difference in the red 
clump levels between the two CMDs. These relative mea¬ 
surements are much more robust than our determination of 
the mean absolute distance modulus above. They are con¬ 
sis tent, at least in a qua litative sense, with the recent work 
of IBalbinot et all (l2015l l. who examined the brightness of 
the red clump in population s lying with in 9° north-west 
of the LMC centre, and also Ivan der Marel k Cionil (l200lh 
who studied the brightness of AGB stars as a function of 
radius and azimuth out to ~ 7°. The distance gradient is 
a result of the LMC geometry, which is usually assumed to 
follow that of an inclined disk with regions to the north and 
east tilted in front of the plane of the sky. 

In the following analysis we aim to use main sequence 
turn-off stars belonging to the dominant old stellar popula¬ 
tion seen in Figure[l]to trace the outer structure of the LMC. 
However, the apparent variation in distance between differ¬ 
ent regions of the galaxy makes uniformly selecting these 
stars problematic because we do not know, a priori , the line- 
of-sight distance to any given point on the sky, and the stel¬ 
lar density is much too low to reliably measure this quantity 
with the necessary spatial resolution. The simplest solution 
is to define a selection box on the combined CMD that is 
sufficiently broad to encompass the full vertical spread in the 
old population. This box is marked in Figure [T] and comes at 
the cost of a slightly higher level of non-LMC contaminants 
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Figure 1 . Hess diagrams for stars in the northern outskirts of the LMC. Left: Stars at all position angles covered by the year one 
DES footprint, and in the radial range 9° — 13°. Centre left: Stars located to the far east of the footprint, at position angles ^ 25°, 
and in the radial range 9° — 13°. Centre right: Stars located to the far west of the footprint, at position angles ^ —55°, and in the 
radial range 9° — 13°. Right: Stars at all position angles covered by the year one DES footprint, and in the radial rang e 7° — 9°. In 
all four panels the magenta line denotes an isochrone from the Dartmouth Stellar Evolution Database ((Potter et al.lf2008r ) with age 11 
Gyr, [Fe/H] = —1.1, and shifted to a distance modulus of 18.35. Note the clear line-of-sight distance variation between the two central 
CMDs relative to this fiducial line. We mark our selection box for old LMC main sequence turn-off stars using a solid line, together with 
a box indicating the location of the red clump belonging to this population. We mark our selection box for exploring contamination from 
the Galactic foreground ("Section 13.3D with a dashed line. In each panel the colour-map represents the stellar density in 0.02 x 0.02 mag 
pixels, smoothed with a Gaussian kernel of a = 0.02 mag, and normalised to 50% of the maximum pixel value. This normalisation has 
the effect of saturating regions of the CMD and introducing artificial fluctuations in the apparent level of non-LMC populations, but 
increases the contrast in low-density regions without moving to a non-linear scale. 


than we might have achieved with a more closely-tailored 
selection regior0. 

We set the faint limit of our selection box at a level that 
mitigates any potential issues due to the onset of (spatially 
variable) detection incompleteness at the level quoted in the 
previous section (see also Section [3731 below"). Assuming that 
the apparent distance gradient from east to west is smooth, 
holding this faint limit constant with position will introduce 
an artificial density gradient on the sky. That is, a patch 
of given stellar density will be observed to be rather less 
dense in the west than in the east due to the CMD selection 
box in the west covering ^0.25 magnitudes less of the old 
population main sequence than it does in the east. We will 
return to this issue later in the paper; for now we simply 
note that the size of this variation is of order ~ 20 per cent 
from the extreme eastern to western edges of the imaged 
LMC region. 


1 Note that the red clump is not included as part of our selection 
region; the box marked in Figure [T] is simply to help facilitate 
our discussion of this feature and its changing level with position. 
Note also that while the red giant branch between the main se¬ 
quence turn-off and the red clump is visible, we do not include 
this region in our selection box in order to minimise contamina¬ 
tion from Galactic foreground stars - see Section 13.31 


3.2 Stellar density map 

In Figure [2] we show a map of the spatial density of stars 
falling in our CMD selection box. This is a tangent-plane 
(gnomonic) projection where the tangent point is located 
at the centre of the LM C (which we assu me to lie a t a = 
82.25°, S — —69.5° as in Ivan der Marel fc Cioni][2001] ). The 
dashed purple lines mark galactocentric radii of 8°, 10°, and 
12°; within the central 8° we have superimposed a wide-field 
optical image of the LMC to help place the periphery into 
perspective. In the DES density map the outer LMC disk 
dominates to a radius of ~ 12° and possesses a mild but 
noticeable non-circular shape, as seen in previous stud ies at 
smaller radial distances (see e.g. Jvan der MarefeOQlh . 

Perhaps the most striking feature on the map is a dis¬ 
tinct stream-like substructure stretching from « 13.5° due 
north of the LMC centre towards the eastern edge of the 
survey footprint. The apparent width of this feature on the 
sky is « 2°. Additional faint structure is visible at commen¬ 
surate radii to the west; based on the density map alone, it is 
not clear whether this constitutes an extension to the main 
body of the substructure or whether it is unrelated (but see 
also Section 13761) . 

The existence of this stream-like feature helps reconcile 
the results of several previous studies of the LMC outskirts: 

(i) A number of authors, while investigating the Carina 
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Figure 2. Map showing the spatial density of stars consistent with being members of the LMC old main sequence turn-off population. 
These were selected from our photometric catalogue using the CMD box marked in Figure [l] The colour-map has been set to saturate at 
a density of 1350 stars deg -2 - this obscures detail in the LMC outer disk, but enhances the visibility of low-density structures without 
moving to a non-linear scale. The bin size is 2' x2 ; ; smoothing has been applied using a Gaussian kernel with a = 3'. On this map, north 
is up and east is to the left as indicated; (£, rf) are coordinates in the tangent plane (gnomonic) projection centred on the LMC. The three 
purple dashed circles mark angular separations of 8°, 10°, and 12° from the LMC centre. Within 8° we display a wide-held optical image 
of the LMC (credit: Yuri Beletsky; used with permission) to help place the size of the periphery into perspective. Our new stream-like 
stellar substructure is clearly evident to the north of the outer LMC disk stretching in the direction of the Carina dwarf, which we mark 
with a magenta point and a surrounding circle that denotes an approximate radius of 1°. The additional faint substructure to the west 
is also indicated. 


dwarf galaxy, have noted the presence of probable LMC stars 
in the foreground. These includ e giants and red clump stars 
observed spe c trosco pically by IMaiewski et al.l |2000) and 
I Munoz et al.l (|2006l ) that appear to constitute a kinemati¬ 
cally cold structure with heliocentric radial velocity ~ 332 
kms -1 , as well as stars at the old main-sequ ence t urn-off 
seen recently in a photometric study bv iMcMonigal et al.l 
(|2014h . The map presented by these latter authors shows the 
LMC populations to be mostly concentrated around ~ 1° 
south of Carina, with a relatively sharp edge. In Figure [2] 
we mark the position of the Carina dwarf with a magenta 
point; the outer circle marks a distance of approximately 1° 
from the dwarf. Extrapolating the position and direction of 
our stream-like feature a short distance off the eastern edge 
of the DES footprint, we posit that these previous studies 
have very likely been detecting members of this substruc¬ 
ture. 

(ii) IMaiewski et al.l (|2009h detected LMC giant stars 
spectroscopically to a radius of ~ 19° to the north-east, 
north, and north-west, and to a ra dius of 22° in the direc¬ 
tion of Carina. On the other hand, ISaha et all ( 2010 ) pho¬ 
tometrically detected main sequence turn-off stars in fields 
to 16° due north of the LMC, but not in fields at 17° and 
19°. Our newly discovered substructure possesses a rather 
sharp drop off in density at ~ 15 — 16° due north of the 
LMC centre. This could explain the apparently conflicting 


results of these two studies if the stellar density beyo nd this 
“edge ” fell belo w the faint det e ction limit of ISaha et all 
(| 201 Oh but not IMaiewski et al.l d 20091 ). LMC populations 
have now been detected photometrically in this direction 
beyon d ^ 15 — 16° (see e.g., Sect ion l3.6l in the present work, 
and iBelokurov fe Koposovl 1 20161 ) . 

3.3 Effects of reddening, incompleteness, and 
contamination 

It is important to check the impact of spatially variable fore¬ 
ground reddening, detection incompleteness, or contamina¬ 
tion of the CMD selection box by non-LMC objects. Each of 
these could, in an extreme case, give rise to an artefact that 
might mimic a stream-like feature; more generally, these is¬ 
sues might alter the observed morphology or properties of 
the substructure. 

The upper left panel in Figure [3] shows the varia¬ 
tion in E(B — V) across the region of interest, from the 
ISchlegel. Finkbeiner fc Davisl ( 1998 ) dust map. As previ¬ 
ously stated, in general the foreground reddening to the 
north of the LMC is mild, with a median value of E(B—V) ~ 
0.05. Although the reddening increases in the east (in the 
direction of the Galactic plane), E(B — V) is still always 
below ~ 0.14 mag. Most importantly, there is no feature in 
the reddening map that is consistent with that of our sub- 
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Figure 3. Maps showing the variation i n the following quantities across the region of interest: Upper left: foreground reddening, derived 
from the lSchlegel. Finkbeiner &; Pavi3 ( 1998) dust map; Upper right: detection incompleteness, calculated in 0.5° X 0.5° boxes; Lower 
left: the spatial density of objects classified as galaxies according to our star-galaxy separation routine; and Lower right: the spatial 
density of stars falling into the foreground selection box defined on the CMDs in Figure [l] For the latter two maps the colour-scaling is 
set to be equal to that in Figure [2] to facilitate easy comparison. In all four panels the three dashed circles denote angular separations of 
8°, 10°, and 12° from the LMC centre, and the Carina dwarf is marked with a magenta point. 


structure. Overall this suggests that variations in foreground 
reddening are not a significant factor for our analysis. 

In the upper right panel in Figure [3] we plot the spa¬ 
tial variation in the r-band magnitude where our catalogue 
starts to become significantly affected by incompleteness. 
We define this in 0.5° x 0.5° bins by locating the turn-over 
in the stellar luminosity function and fitting a curve to de¬ 
termine the level at which the number counts fall to 75% 
of this peak value. The intrinsic bin-to-bin scatter in this 
calculation is roughly ±0.1 mag. In general the complete¬ 
ness level is relatively consistent across most of the region 
of interest, although there is a systematic move to brighter 
magnitudes in the east. This could be due to a general de¬ 
crease in image quality or the effective integration time (i.e., 
fewer available epochs) near the edge of the DES year one 
survey footprint. There is also a noticeable move to brighter 
magnitudes inside 8° of the LMC centre due to crowding. 
Near Carina there is a small patch where the completeness 
level moves to much brighter magnitudes - this is a region of 
poor data quality (see below). Importantly, apart from this 
small patch, variations in the completeness level always oc¬ 
cur at substantially fainter magnitudes than the cut-off level 
in our CMD selection box. This is by design, to ensure that 
spatially variable incompleteness does not affect our analy¬ 
sis. The bad patch of data falls squarely on the eastern end 


of our substructure, and may explain why, in Figure [2] the 
feature appears to peter out towards the edge of the survey 
area. 

Despite our CMD selection box being optimised to iden¬ 
tify LMC main sequence turn-off stars, there will inevitably 
be some level of contamination by non-LMC objects. The 
two primary sources of contamination are unresolved back¬ 
ground galaxies, and stars belonging to the Galactic fore¬ 
ground. To investigate the former, we constructed a density 
map of objects classified as galaxies rather than stars ac¬ 
cording to the star -galaxy separat i on ro utine described in 
Section [2] (see also IKoposov et alJ 12015) . While these ob¬ 
jects are of course not unresolved galaxies, the expectation 
is that they will behave in a similar fashion. This density 
map is shown in the lower left panel in Figure [3] In general, 
the distribution of galaxies is relatively uniform across the 
survey area, with a few filamentary regions of excess. The 
substantial excess towards the central regions of the LMC 
is due to crowding and the consequent misclassification of 
stars as galaxies. The patch of bad data near Carina reap¬ 
pears in this map - poor image quality in this region results 
in both the misclassification of stars as galaxies, and the 
previously-observed bright completeness limit. 

Although there are a variety of Galactic stellar popu¬ 
lations in the foreground, not all of these will contaminate 
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our LMC selection box on the CMD. In Figure [I] we mark a 
second selection box (using a dashed line), that is designed 
to pick out old main sequence turn-off stars in the magni¬ 
tude range 21 < r < 19 - i.e., at distances of ~ 10 — 30 kpc. 
The main sequences associated with these turn-off popula¬ 
tions cut back towards the red and across the upper part of 
our LMC selection box; therefore, they serve as a good proxy 
for exploring how the foreground contamination should vary 
with position. We show a density map for stars falling in this 
secondary selection box in the lower right panel in Figure [3] 
Unsurprisingly, since most of the stars in the distance range 
under consideration will belong to the Galactic halo, the dis¬ 
tribution is very smooth. There is only a very mild gradient 
towards the east, approaching the Galactic plane, presum¬ 
ably due to the nearest contaminants. In general the ex¬ 
pected level of contamination from foreground populations 
is only a small fraction of that of LMC stars, even in the 
outskirts of the LMC disk and including our new stellar 
substructure. We note that this analysis does not include 
main sequence turn-off stars at distances beyond ~ 30 kpc, 
which will also fall into the LMC selection box; however in 
these regions the Galactic halo is expected to have a negli¬ 
gible surface density even compared to the foreground map 
plotted in Figure [3] Finally, we note an interesting apparent 
overdensity of foreground stars towards the centre of the 
LMC - this excess is caused by young LMC populations, 
that are present only within « 8 — 9° of the centre, cutting 
across the foreground selection box on the CMD (see the 
rightmost panel in Figure [lj. 

In summary, we have shown that spatially variable fore¬ 
ground reddening, detection incompleteness, and contami¬ 
nation of the CMD selection box by non-LMC objects can¬ 
not explain the presence of the stream-like feature observed 
in Figure [21 and it is thus reasonable to conclude that this 
is a genuine low surface-brightness substructure in the out¬ 
skirts of the LMC. 

3.4 Properties of the stream-like feature 

Figure [4| shows a second version of our stellar density map, 
tuned to enhance low surface brightness features in the LMC 
periphery. The faint structure in the outskirts of the disk to 
the west is clearly visible, as is some degree of structure 
to the south of the main stream-like feature stretching to 
the east. In order to explore the properties of the stellar 
populations belonging to this feature, we define a selection 
box around its densest section. As marked in Figure [U this 
box is split into two equal regions - one covering the eastern 
half of the substructure and one covering the western half. 

In the left-hand panel of Figure [5] we display a Hess 
diagram showing all stars falling within the selection box 
(i.e., within both halves). Although the foreground contam¬ 
ination is relatively severe, the signature of the substructure 
is clearly evident on this CMD - both an old main sequence 
turn-off and a red clump are visible. The fiducial isochrone 
that we used in Figure [I] appears to match the locus of sub¬ 
structure stars quite closely, suggesting that the stellar pop¬ 
ulation (s) that constitute this feature are very similar to the 
dominant old stellar population seen in the outer LMC disk. 

To investigate this further, we define a second box im¬ 
mediately to the north of the selection box for the substruc¬ 
ture, that covers an equivalent area of sky and which we use 



£ (degrees) 


Figure 4. A second realisation of our stellar density map, with 
the colour-scale set to saturate at 900 stars deg -2 and an in¬ 
creased degree of smoothing (using a Gaussian kernel of width 
a = 8 7 ). The selection box used to explore the properties of the 
stellar populations in the stream-like substructure is marked with 
a solid line, and the corresponding region used for the field sub¬ 
traction is marked with a dashed line. As before, the dashed circles 
indicate angular separations of 8°, 10°, and 12° from the LMC 
centre. 

to define a Hess diagram for the foreground populations that 
contaminate our stream CME0. We smooth both the “sub¬ 
structure” and “field” Hess diagrams with a Gaussian kernel 
of a — 0.05 mag and then subtract the field distribution from 
that for the substructure via a straight grid point by grid 
point differencing. The resulting field-subtracted Hess dia¬ 
gram is displayed in the centre-left panel in Figure [5j The 
fiducial isochrone clearly provides an excellent description 
of the substructure populations; indeed even the distance of 
these populations appears to match quite closely to that of 
the isochrone. The red clump is clearly visible and corrobo¬ 
rates this conclusion. 

We repeated this procedure for stars in the western and 
eastern halves of the substructure selection box, to obtain 
the Hess differences plotted in the centre-right and right- 
hand panels of Figure [5] No significant variations in distance 
are evident from west to east along the feature, although we 
note that the CMD for the eastern half is very noisy due 
to the increasing level of foreground contamination in this 
di rection. This resu l t is co nsistent with the measurements 
of lMcMonigal et"ahl (|2014l b who found that the LMC stars 
they observed in front of the Carina dwarf lie ~ 46 kpc away 
(i.e., at a distance modulus ~ 18.3). 

The arc length along the substructure from 13.5° due 
north of the LMC centre to 1° due south of the Carina dwarf 
is ~ 12.5° on the sky, corresponding to a physical length of 
~ 10 kpc. To estimate a lower limit for its luminosity, we 
added up the total amount of light from stars in both parts 
of our selection box that also sit within the main sequence 
turn-off selection box on the CMD, and subtracted the con¬ 
tribution of non-stream members using the field selection 
box. Ignoring any differences between DES g , r and SDSS 

2 Note that is is likely that this box contains some LMC-like stars 
(see Sect ion [3Tl) : however, these have a much lower spatial density 
than the substructure populations and so are of no consequence 
here. 
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Figure 5. Hess diagrams illustrating the properties of stars in the stream-like overdensity. Far left: Stars lying in both halves of the 
selection box marked on the map in Figure|4] Centre left: Field subtracted version of this Hess diagram. Centre right: Field subtracted 
Hess diagram for stars lying in the western half of the selection box. Far right: Field subtracted Hess diagram for stars lying in the 
eastern half of the selection box. In the left-most panel the CMD is binned into 0.02 x 0.02 mag pixels and smoothed with a Gaussian 
kernel of a = 0.02 mag. In the other three panels the CMDs are binned into 0.04 x 0.04 mag pixels and smoothed with a Gaussian kernel 
of a — 0.05 mag. In all four panels the colour-scale is set to saturate at 40% of the maximum pixel value, and the same isochrone fiducial 
from Figure ^ is marked with a magenta line. The same main sequence turn-off and red clump selection boxes as in Figure H are a l so 
plotted. 


g , r, which are known to be small (e.g.. iTucker et aI1ll2007 h 
we transformed from SDSS g and r magnitudes to the John¬ 
son E-band using the relationship defined by Lupton on the 
SDSS DR12 web pagj^to obtain Mv = —5.7. 

However, the main sequence turn-off stars sitting in our 
CMD selection box represent only a fraction of the total sub¬ 
structure population - for example red giants, red clump 
stars, and faint main sequence stars are all excluded. To de¬ 
termine an approximate correction to account for light due 
to these missing populations we generated a luminosity func¬ 
tion from the Darthmouth Stellar Evolution Database with 
the same age and meta llicity as our fiducial isochrone, and 
based on an underlying iKrou toI (l200lh mass function. This 
revealed that our CMD selection box encompasses only « 20 
per cent of the total light expected from a stellar popula¬ 
tion with these characteristics, implying a total luminosity 
of My ~ —7.4. This is still, of course, a lower limit to the 
true luminosity of the structure because our spatial selec¬ 
tion box covers only its brightest portion. We also note pos¬ 
sible systematic uncertainties of up to ~ 0.5 mag due to any 
spread in the characteristics of the substructure populations 
that we are unable to detect using the presently-available 
CMDs, as well as the inclusion of underlying disk popula¬ 
tions that may or may not be appropriate depending on the 
origin of the feature (see Sections 13.61 and l4l). Assuming that 
My ~ —7.4 , the mean apparent E-band surface brightness 


of the substructure within the selection box is ~ 31.8 mag 
arcsec -2 . 


3.5 Structure of the outer LMC 

In order to try and obtain some insight into the origin of the 
stream-like overdensity, we now turn to a brief investigation 
of the structure o f the outskirts of the LMC. Following previ¬ 
ous s tudies fe.g. JWeinberg &; Nikolaevll200ll : [Balbinot et al.l 
l2015l i we first attempt to model the geometry of this region 
using a circular exponential disk that is in clined to the p lane 
of the sky. We adopt the formalism of lvan der Marel &; Cionil 
(l200lh and define coordinates (p, (j)) on the celestial sphere 
for a given point with right ascension a and declination J, 
such that 

cos p — cos J cos Jo cos (a; — ao) + sin J sin Jo 


tan 0 — 


cos J sin Jo cos(a — ao) — sin J cos Jo 
cos J sin(a — ao) 


(i) 


where (ao, Jo) are the right ascension and declination of 
the centre of the LMC. The coordinate p is the angular dis¬ 
tance between (a, J) and (ao, Jo), while 4> is the position an¬ 
gle of (a, J) with respect to (ao, Jo) taken counterclockwise 
from the west@. As described earlier, we assume ao = 82.25° 
and Jo = —69.5°. This represents the approximate central 


3 http://www.sdss.org/drl2/algorithms/sdssUBVRITransform/ 


4 This is different from the usual astronomical convention where 
position angle is measured counterclockwise from the north. 
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point of the outer n umber d ensity contours measured by 
Ivan der Marell (|200lh (see also Ivan der Marel k, Cionill2001f ) 
and as such is appropriate for the present problem. 

We next introduce a Cartesian coordinate system (x, y ) 
in the plane of the LMC dislQ; for simplicity we assume that 
the disk is infinitely thin. The disk plane is inclined with 
respect to the sky by an angle i around an axis (the line of 
nodes) that sits at position angle 0 , which is again measured 
counterclockwise from the west. The line of nodes represents 
the intersection of the LMC disk plane with the plane of the 
sky. The coordinates (x, y ) may be calculated from (p, (f>) 
according to 

Do cos i sin p cos(0 — 0) 
cos i cos p — sin i sin p sin(0 — 0) 

_ _ Do sinpsin(0 - 6) _ 

^ cos i cos p — sin i sin p sin(0 — 0) ’ 

where Do is the distance to the centre of the LMC, which, as 
noted previously, we assume to be 49.9 kpc. The line-of-sight 
distance D to any point on the LMC disk plane is given by 


cos i cos p — sin i sin p sin(0 — 0) 

In our model the surface density in the plane of the 
LMC disk varies as po exp(— r/r s ), where the radial distance 
r 2 = x 2 + y 2 , the scale length is r s , and the central surface 
density is po- Note that we do not include a background 
contribution as the density of contaminants is much lower 
than the density of LMC stars over the region where we fit 
our model (see below, and Section 13.61) . We divide the tan¬ 
gent plane (£, 77 ) into a grid of 10 / x 10 / cells and compare 
the observed number of stars per cell to that predicted by 
the model. In doing so, there are three additional factors we 
need to take into account. First, the gnomonic projection 
(a, (5) —»• (£, p) does not conserve area. The distortion is ap¬ 
proximately proportional to l/(cos p) 3 , such that a constant 
cell size on the tangent plane corresponds to an increas¬ 
ingly smaller area on the celestial sphere as the separation p 
from (a 0 , £ 0 ) becomes larger; we use this to correct the pre¬ 
dicted star counts. Second, we must account for the fact that 
the line-of-sight distance to the inclined LMC plane varies 
across the sky. A cell for which the distance is shorter spans 
a smaller physical area on the plane than does a cell for 
which the distance is larger; the necessary scaling of counts 
per cell thus varies as D 2 . 

Finally, as previously described in Section 13.11 the line- 
of-sight distance variation also tends to move stars in and 
out of our fixed CMD selection box, which can lead to arti¬ 
ficial density gradients across the sky. To attempt to incor¬ 
porate this into our model we define an empirical correction 
function A(0), as follows. We take all stars in the radial range 
9° — 13°, as plotted in the left-hand panel of Figure [2 and 
split them by position angle into six equal-sized subgroup^] 
(the central-left and central-right panels in Figure [I] repre¬ 
sent the extreme eastern and western of these ensembles). 
For each subgroup we determine the vertical shift necessary 

5 Note that this is th e (x', y') coordinate system defined by 
Ivan der Marel Sz Cionil (|20011 ). 

6 This represents the finest partitioning for which we maintain 

adequate statistics on the Hess diagrams. 


to align the fiducial isochrone with the ridgeline on the Hess 
diagram, relative to the shift required for the eastern-most 
segment. We also determine the relative offsets in the red 
clump levels by measuring the position of the peak of a his¬ 
togram in r magnitude within the red clump selection box, 
and find strong consistency with the offsets derived using the 
fiducial isochrone. We reiterate that these relative distance 
measurements are much more robust to potential system¬ 
atic uncertainties than absolute distances determined from 
isochrone fitting. 

The upper panel of Figure [6] shows the results of our 
measurements (recall that the position angle 0 is measured 
counterclockwise from the west). The line-of-sight distance 
appears to increase smoothly and monotonically from east to 
west, and the gradient becomes shallower in the far west; the 
ove rall shape of the c u rve is consi stent with being sinusoidal 
(cf. Ivan der Marel & Cionil l200lh . Because our data cover 
only about one-third of the outer LMC disk we are not able 
to constrain the phase and amplitude of the sinusoid very 
precisely; however we overplot, relative to the distance zero- 
point at 0 « 125°, a representative fit which has a peak-to- 
peak amplitude of 0.34 mag and a phase shift of 275°. Note 
that the stars in each azimuthal bin marked on the plot have 
a different mean angular distance from the LMC centre, due 
to the elliptical shape of the LMC disk on the sky. 

To determine how the distance gradient affects our star 
counts, we took all stars found inside the CMD selection 
box for the eastern-most group, increased their magnitudes 
by the offset appropriate for each of the other segments in 
turn, and then counted how many still fell within the CMD 
selection box. The outcome is shown in the lower panel of 
Figure [6j As expected there is a substantial decrease in the 
number of stars falling in the selection box as a function of 
position angle. The magnitude of the variation is of order 20 
per cent from east to west across the survey footprint. Our 
measurements are well fit by a second-order polynomial of 
the form 

H&) = 0.1255 (yqo) + 0.3418 (^) + 0.76 . (4) 

Returning to our inclined exponential disk model, the 
overall number of stars predicted in a given cell j is 

mj(p o, r s , 0, i) = p 0 D 2 (cos pj) 5 3 X^j) exp(-r j /r s ). (5) 

The probability of finding the observed number of stars rij 
in this cell is given by a Poisson distribution 

m nj e - rn 3 

p A n i \ m j) = —— j - , (6) 

lljl 

so that the logarithm of the likelihood of the data across all 
cells is 

£{n, p, <j> \no, r s , 0, i) = InH-Pl 

3 

= Y ln ( m i) _ m i ~ ln (% ! )] • ( 7 ) 

3 

7 Although the locus of points appears to share a very similar 
shape to that for the relative distances in the panel above, such 
that the functional form may well be formally sinusoidal, a poly¬ 
nomial fit works well over the range in azimuth considered here. 
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0 (deg) 

Figure 6. Upper panel: Relative distance measurements as a 
function of position angle cf) for the six segments of the outer LMC 
disk (radii 9° — 13°) described in the text. Recall that cf) is mea¬ 
sured counterclockwise from the west. All distances are measured 
relative to the easternmost segment. Black points mark distances 
derived from aligning the fiducial isochrone with the CMD ridge¬ 
line in the main sequence turn-off region, while red points mark 
distances derived from the red clump. There is good agreement 
between the two sets of measurements. The black solid line shows 
a sinusoid with a peak-to-peak amplitude of 0.34 mag and a phase 
shift of 275°. Lower panel: Star counts for each of the six seg¬ 
ments relative to the easternmost segment. The variation is of or¬ 
der 20% from east to west; the measurements are well described 
with a second-order polynomial (marked with the red line). 


The model includes implicit information about the distance 
to the stars in each cell, as this affects the observed den¬ 
sity through both the D 2 and A(0) terms. Ideally we would 
also self-consistently include an explicit comparison of the 
measured distance to a given cell with that predicted by the 
model; however, as previously noted the stellar density is too 
low to allow distance measurements with useful precision at 
the necessary spatial resolution. 

We use the M arkov chain Monte Carlo sampler emcee 
(|Foreman-Mackev et al.ll2013l ) to infer the most likely set of 
model parameters given our observations and the assump¬ 
tion of flat priors. We restrict /io and r s to be positive, and 
0° ^ i < 90°. With this bound on the inclination, and our 
observation that LMC stars lying to the east of the survey 
footprint are closer than those lying to the west, it follows 
that 180° ^ 0 < 360°. Our calculation considers all cells 
with 9° ^ p ^ 13°, where the inner bound is set to min¬ 
imise contamination in our CMD selection box from young 
populations (which may not b e well descr ibed by an ex¬ 
ponential disk model - e.g., iBalbinot et al.ll2015h . and the 
outer bound ensures we do not encroach on the region occu¬ 
pied by the substructure (which is a clear deviation from an 
exponential disk). We use an empirically-determined mask 
to excise the ~ 1% of cells that intersect with small gaps in 


the spatial coverage, regions of bad data, or the two LMC 
globular clusters that fall in the survey footprint (Reticu¬ 
lum and ESO121-SC03). Figure [7] (upper panel) shows the 
one- and two-dimensional marginalized posterior probabil¬ 
ity distributions for our model parameters, and Table [I] lists 
the most-likely values for these quantities. For consistency 
with the literature, in which the position angle of the line of 
nodes is usually taken co unte rclockwise from the north, we 
follow Ivan der Marel &; Cionil (|200lh and report our results 
in terms of © = 0 — 90°. The formal uncertainties for all 
parameters are very small, showing that the algorithm has 
converged well; however it is likely that the true uncertain¬ 
ties, including systematics, are substantially larger. 

Our inferred values for the inclination and scale length 
are consistent with estimates from previous work on the 
outer LMC disk, although some degree of variation is seen 
between measurements in the literature due to the different 
regions covered by various studies, and t he use of differ- 
ent tr acer populations. For exam ple Jvan der Marel &; Cionil 
(|200lh and Ivan der Marell (|200lh find i = 34.7° ± 6.2° and 
r s « 1.3 — 1.5 kpc by tracing AGB and R GB stars ov er 
the radial range ~ 2° — 7°, while I Weinberg fc Nikolaevl 
(l200lh find i = 24.0° d= 0.3° and r s = 1.42 ± 0.01 kpc 
over a similar radial ra nge usin g a variety of old and inter¬ 
mediate age tracers. I Saha et ah ( 2010 1 measure r s — 1.15 
kpc by tracing old main sequence turn-off stars in a stripe 
spanning ~ 6° — 16° d ue north of the LMC a nd assum¬ 
ing the geometry o f van der Marel &; Cionil (|200lh . More re¬ 
cently, IBalbinot et al ] (l2015h obtain i = 32.9° ± 0.4° and 
r s = 1.41 d= 0.01 kpc by tracing old main sequence turn¬ 
off stars in a region spanning a similar radial range to the 
north-east. 

The lower panel of Figure 0 shows our model plotted on 
the LMC density map. There is good general correspondence 
between the shape of the model and the observed isodensity 
contours. Our inferred value for O implies the line of nodes 
lies ~ 20° from vertical on the map. In Figure [8] we plot the 
line-of-sight distance profile for the best-fit model (blue line). 
Although the model does exhibit a clear distance gradient 
from east to west, the gradient is too steep in the sense that 
the magnitude offset between the far east and far west of 
the survey area appears too large. 

There are two reasons for this discrepancy. First, the 
position angle of the line of nodes is not a good match for the 
observed distance profile. The sinusoid plotted in Figure [6] 
implies that the line of nodes should sit at an angle ~ 5° east 
of north; however that for our inclined circular disk model 
is oriented 25° further westwards. Forcing the line of nodes 
in the model to match that implied by the distance profile 
would move the blue curve left and downwards on the plot 
(recall that all our distance measurements are relative to the 
easternmost segment). However, this would then mean that 
the shape of the model would not fit the isodensity contours 
so well - for an inclined circular disk, the orientation of the 
line of nodes must coincide with the direction of the major 
axis of the (elliptical) isopleths on the sky. 

Second, regardless of the orientation of the line of nodes, 
the extremes of the line-of-sight distance in our model (at 
any azimuth, not just within the survey footprint) over an¬ 
gular separations from the galaxy centre commensurate with 
those for our six segments, span a range of ~ 0.45 mag in 
distance modulus, which is ~ 30 per cent larger than the 
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Table 1. Most-likely parameter values for our two inclined exponential disk models. 


Model 

Radial range 

yo 

r s 

© 

i 

*p e 



(stars cell -1 kpc -2 ) 

(kpc) 

(deg) 

(deg) 

(deg) 

Circular disk 

9° - 13° 

20.95 ± 0.26 

1.50 ±0.01 

159.59 ±0.12 

33.14 ±0.09 


Elliptical disk 

9° - 13° 

12.80 ±0.25 

1.64 ±0.01 

(185.0) 

25.18 ±0.71 

34 . 68 t 2.24 °' 092 ± 0-001 


peak-to-peak amplitude suggested by our observed distance 
profile. This implies that the best-fit inclination is too high; 
however, a reduction in the inclination would again mean 
that the shape of the model on the sky would not fit the 
isodensity contours so well, as this parameter controls the 
degree of ellipticity on the sky. 

Retaining the assumption of a planar geometry, one way 
of resolving these issues is if the outer LMC disk is intrin¬ 
sically elliptical rather than circular, perhaps d ue to tida l 
stresses acting on the disk (see e.g., Ivan der Marell [200lh . 
To test the feasibility of this hypothesis we replaced the cir¬ 
cular disk in our model with an ellipse defined on the (x, y) 
plane by 

(x cos'ip — y sin ip) 2 (xsinvp ± ycos'ip) 2 , . 

-2- 1 -oTi-\2- = 1 • (°) 

a z a z { 1 — e) z 

The ellipse has semi-major axis a and ellipticity e— (1 — q) 
where q is the axis ratio, and is rotated on the plane by 
an angle 'ip which, for consistency with <p and 6 , we mea¬ 
sure counterclockwise from the “west” (i.e., when 0 and i 
are zero, positive 'ip gives rotation north of west on the sky). 
The scale length is defined along the major axis of the el¬ 
lipse on the disk plane. The extra parameters in this model 
mean that the orientation of the density contours on the sky 
is decoupled from the position angle of the line of nodes; 
moreover, the apparent ellipticity of the isopleths is at least 
partly controlled by the intrinsic shape of the elliptical disk 
rather than solely by the inclination. 

We tried to fit the full six-parameter elliptical model 
to our data; however because the survey footprint covers 
less than 40 per cent of the outer LMC disk, several of the 
parameters are largely degenerate and the algorithm fails to 
converge. We thus attempt to simplify the problem by fixing 
the orientation of the line of nodes to the value 0 = 185° 
implied by our distance profile - essentially we force the 
phase of the sinusoidal distance variation to agree with the 
observations and then search for a solution that matches the 
shape of the density isophotes on the sky. 

The upper panel of Figure [9] shows the marginal¬ 
ized posterior distributions for this five-parameter prob¬ 
lem; the most-likely values of the parameters are listed 
in Table |T] The scale length is similar to that for our 
circular model. The inclination is substantially lower, al¬ 
though it is still well within the range of inclination mea¬ 
surement s_JoundJn_t he literat ure at smaller radii (see e.g., 
ISubramanian fe Subramaniaml l201(j and the compilation 
therein) and at larger radii as discussed above. Interestingly, 
our value for i matches quite closely with a linear extrapola¬ 
tion of the radial variatio n in this quantity seen out to ~ 7° 
by Ivan der Marel &; Cionil (l200lh (their Figure 8). Finally, 
the intrinsic ellipticity of our model disk is somewhat lower 
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Figure 7. Upper panels: Marginalized one- and two- 
dimensional posterior probability distributions for the four pa¬ 
rameters in our inclined circular exponential disk model. These 
are the central density yo, the scale length r s , the position an¬ 
gle of the line of nodes © (measured counterclockwise from the 
north), and the inclination angle i. The contours show the lcr, 2cr, 
and 3cr confidence levels. Lower panel: The best-fitting model 
overplotted on our stellar density maps. Various isodensity con¬ 
tours are marked. The model is plotted at four different radii in 
the disk plane - 8, 9, 10, and 11 kpc. The colour-scaling indicates 
the line-of-sight distance to each point on the model, as indicated. 
In general there is good correspondence between the shape of the 
model and the isopleths. North is up and east to the left. 
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0 (deg) 

Figure 8. Relative distance measurements as a function of po¬ 
sition angle 0 f° r our six outer disk segments, along with the 
profiles for our inclined circular and elliptical disk models (blue 
and magenta lines, respectively). The black points represent the 
mean value of the main sequence turn-off and red clump distance 
measurements for each segment. 


than that of Ivan der Marell ( 2001 ). which has e « 0.31 out 
to - 7°; indeed our model is not far from being circular. 

The bottom panel of Figure [9] shows our elliptical disk 
model plotted on the LMC density map. As with our circular 
model, there is in general very good agreement between the 
shape of the model and the observed isodensity contours. 
The line-of-sight distance profile for the elliptical model is 
plotted with a magenta line in Figure [ 8 ] This is an excellent 
match to the observed distances. In particular, the extremes 
of the line-of-sight distance in the model at galactocentric 
radii comparable to those for our six segments span a range 
of ~ 0.34 mag in distance modulus, which is entirely consis¬ 
tent with the peak-to-peak amplitude inferred observation- 
ally. Although we fixed the position of the line of nodes to 
match the observations, it is encouraging that the correct 
amplitude of variation arises naturally from the maximum 
likelihood algorithm even though this is simply trying to 
match the observed stellar density on the sky without in¬ 
corporating any explicit distance dependence. 

In summary, while an inclined circular disk model can 
describe the observed density distribution of the LMC out¬ 
skirts on the sky, it cannot simultaneously reproduce the 
observed azimuthal line-of-sight distance profile. To do this 
successfully requires a more general inclined elliptical disk 
model. Even so, there is still one prominent discrepancy be¬ 
tween this model and our observations. Equation [3] suggests 
that the line-of-sight distance to points on the line of nodes 
should be D — Do /cos p « 18.52, whereas we measure a dis¬ 
tance modulus of ~ 18.35 in the relevant azimuthal segment. 
While there are certainly systematic effects to be concerned 
about in the derivation of this distance, these would have 
to be quite substantial (^0.15 — 0.2 mag) in order to com¬ 
pletely account for this offset. One possibility is that the 
LMC geometry becomes non-planar at large galactocentric 
separations, and indeed a number of previous work s have 
seen evidence for such deviations (e. g. Nikolaev et al.ll2 004: 
ISubramanian &; Subramaniaml | 20 10 ; Ba lbinot et al.l l2015l ) . 


Flaring and/or warping of the outer LMC disk could ex¬ 
plain the difference between our absolute distance scale and 
that predicted by the elliptical disk model. 

Twisting of the line of nodes with galactocentric sep¬ 
aration is a key signature expected if the LMC disk is 



Figure 9. Same as Figure [7| except for our inclined elliptical ex¬ 
ponential disk model. The five parameters are the central density 
po , the scale length along the major axis r s , the inclination an¬ 
gle z, the position angle of the elliptical disk on the LMC plane 
0, and the ellipticity of the disk e. In the lower panel the model 
is plotted at radii of 8.5, 9.5, 10.5, and 11.5 kpc in the plane of 
the disk. Again, there is good general correspodence between the 
shape of the model and the isodensity contours. As previously, 
north is up and east to the left. 


warped (e.g., iNikolaev et al.ll2004l ). Most studies of the in¬ 
ner LMC infer a position angle for the line of nodes to be 
in the range ~ 140° — 170 ° east o f north (s ee e.g ., the com¬ 
pilation in Subramani an &; Subramaniaml I2OI0I ). although 
Ivan der Marel fe Cionil (| 200 ll ) found a value of 122.5° d= 8.3° 
east of north out to ~ 7°. Our distance profile suggests 
© ~ 185° at large galactocentric radii, and is thus perhaps 
indicative of just such a twist. 

It is also worth noting that the outermost density con¬ 
tours in the lower panels of Figures [7] and [9] deviate notice¬ 
ably from a purely elliptical shape. This does not just occur 
in the region of the stream-like overdensity - for example, 
there is a clear outwards bulge in the contours to the north¬ 
west, near (£, 77 ) = (—10, 7). In the left panel of Figure flOl 
we show a logarithmic stellar density map of the outer disk, 
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Figure 10. Spatial density maps of the outer LMC disk. Left panel: A map similar to that shown in Figure [2] but with a logarithmic 
scaling to reveal features in the region interior to ~ 13°. The LMC globular clusters Reticulum and ESO121-SC03 are clearly visible, 
at (£, rf) « (—7, 10), and « (4, 8), respectively. Also apparent are several possible deviations from a smooth distribution, primarily in 
the form of sharp edges near (£, rf) « (—10, 7) and « (6, 10). The two dashed circles indicate radii of 9° and 13° from the LMC centre, 
corresponding to the range over which our planar disk models were fit. Right panel: A map showing, with a linear scale, residual star 
counts once our best fitting elliptical disk model has been subtracted, together with a uniform background (see Section [3761 below). In 
addition to our remote stream-like feature, widespread deviations from a smooth planar disk are clearly evident between 9° — 13°. Within 
9° the model does not provide a good description of the data due to the presence of additional young stellar populations in our CMD 
selection box. 


revealing a sharp edge in this region, as well as additional 
edges at other azimuthal angles and radii - for example, near 
(£, rf) = (6, 10). These features suggest that not only is the 
outer LMC disk possibly distorted, its structure may not be 
smooth. 

To explore this idea in more detail, we construct in the 
right hand panel of Figure [10] a map showing the residual dif¬ 
ferences between the observed stellar density and our best¬ 
fitting smooth elliptical disk model. In the region over which 
the best-fitting model was calculated (i.e., at radii between 
9° — 13° from the LMC centre), substantial positive residu¬ 
als are evident. These take the form of a low suface bright¬ 
ness ring-like structure spanning all azimuthal angles in the 
survey footprint, but densest to the west. The sharp-edged 
features seen in the logarithmically scaled map of the disk 
are due to this structure. Its nature is not clear from the 
presently-available data - the structure could represent an 
overdensity in the disk itself, akin to a spiral arm, or, more 
speculatively, it may represent a stellar stream of some de¬ 
scription, perhaps due to an accreted satellit^E Further in¬ 
vestigation is clearly warranted. 

Inside 9° our model does not provide a good fit to the 
data. This is largely due to the presence of young stellar 
populations at such radii, as discussed in Section [3.11 Out¬ 
side 13° the model does an excellent job of describing the 
structure of the LMC disk, at least to the east of £ = 0, even 
though data at such large radii were not used to calculate 
the model. To the west of £ = 0 our new stream-like feature 
is clearly visible stretching towards the Carina dwarf. The 
map clearly shows that this feature “terminates” due north 
of the LMC- further east and further south the data are well 
described by our planar disk model. 


8 It is interesting that the peculiar remote LMC globular clu sters 
Reticulum and ESO121-SC03 (e.g., iMackev et ahll2004l . [20061) sit 
on either edge of the structure. 


3.6 Outer surface density profile and the extent 
of the LMC 

In the upper panel of Figure [lT] we plot the radial surface 
density profile for stars in our main sequence turn-off CMD 
selection box, where the radius is measured along the major 
axis in the plane of the elliptical disk from the best-fitting 
model described above. We corrected all star counts using 
the empirical relation in Equation [4] The predictions of the 
model (including a background contribution, see below) are 
marked with a red line. This provides an excellent fit over 
the radial range we used for the input data (~ 9° — 13°). 
This is despite the presence of the low surface brightness 
structure identified in Figure [TO] which is disguised by the 
azimuthal averaging inherent in the construction of Figure 
rm plus the logarithmic scale. 

As previously noted, for radii inside ~ 8.5 kpc (« 9°) 
the density is higher than the model, due to increasing con¬ 
tamination from younger LMC populations. Outside ~ 12.5 
kpc (~ 13.5°), the density again becomes noticeably higher 
than predicted by the model, due to the presence of the 
stream-like overdensity and contamination of the CMD se¬ 
lection box (primarily by foreground sources). At very large 
radii, the density profile asymptotes to an approximately 
constant background level (marked with a horizontal dashed 
line); we discuss this in more detail below. 

The lower panel of Figure [TT] shows the radial surface 
density profile split into two - blue for stars sitting to the 
east of a line running due north from the LMC centre (i.e., all 
stars east of £ = 0), and magenta for stars sitting to the west 
of this line. Although the profiles match quite well at inner 
radii, they deviate substantially outside ~ 12.5 kpc (i.e., 
~ 13.5°). We mark this radius with a vertical dotted line in 
the plot. This deviation is the signature of the main body of 
our new stream-like overdensity, which lies almost entirely 
to the east of £ = 0. Outside 12.5 kpc, the density profiles 
sit apart for almost the full extent of our measurements. 
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Figure 11. Radial surface density profiles for stars falling in our 
main sequence turn-off selection box on the CMD. Upper panel: 
Density profile for all stars. Lower panel: Density profiles for 
stars lying to the east and west of £ = 0 (blue and magenta 
points, respectively). In both panels the prediction of our best¬ 
fitting elliptical disk model, plus the background, is marked (red 
line in the upper panel; black line in the lower panel). Also marked 
in both panels is the background level (black dashed line) as de¬ 
termined empirically using the western profile, and the four key 
radii on the LMC disk plane discussed in the text (vertical dotted 
lines). These are: (i) 12.5 kpc (w 13.5°), which marks the inner¬ 
most edge of the main stream-like overdensity and the radius at 
which the eastern and western profiles begin to diverge; (ii) 14.5 
kpc (« 15.7°), which marks the outer edge of the faint features 
in the region to the west of £ = 0; (iii) 16.5 kpc (« 17.8°), which 
marks the outer edge of the main stream-like overdensity in the 
region to the east of £ = 0; and (iv) 18.5 kpc (« 20°), which 
marks the radius beyond which we do not see evidence for LMC 
populations in any direction. 


We interpret these profiles with the aid of our stellar 
density map, which is presented in Figure fl2l with a high de¬ 
gree of smoothing and a colour scaling that emphasises very 
low-density features at large radii. The innermost dashed el¬ 
lipse on this map corresponds to a radius of 12.5 kpc; it is 
easy to see why the profiles deviate so strongly immediately 
outside this bound. The western profile (magenta) decreases 
smoothly out to a radius of ~ 18.5 kpc, where it flattens to a 
constant background. This radius is marked with the outer¬ 
most dashed ellipse on the density map, and the outermost 
vertical dotted line on the profile. The density map shows 
that to the west of £ = 0 and between 12.5 kpc and 18.5 kpc, 
there are clearly still LMC stars present. Most of these re¬ 
mote stars sit within ~ 14.5 kpc (also marked), which traces 
the approximate outer edge of the faint western structure 
that we remarked upon in Section ET21 however, there are 
also even fainter features present in the range 14.5 — 18.5 
kpc. These extremely low surface density populations are 
likely to be associated with the LMC as they broadly fol¬ 
low the same elliptical shape on the sky. As noted in our 



£ (degrees) 


Figure 12. A final incarnation of our stellar density map, with 
the colour-scale set to saturate at 600 stars deg -2 and a large 
degree of smoothing (using a Gaussian kernel of width cr = 10'). 
As before, north is up and east to the left. The dashed ellipses 
correspond to the four key radii plotted in Figure fTTI and described 
in the text. The innermost ellipse marks a radius in the disk plane 
of 12.5 kpc, while the outermost ellipse marks 18.5 kpc. Between 
these inner and outer ellipses, a radius of 14.5 kpc is marked in 
grey to the west of £ = 0, and a radius of 16.5 kpc is marked in 
grey to the east. For convenience the line £ = 0 is also plotted. 


discussion of Figure ITOl above, in this region to the west of 
£ = 0 our exponential disk model, which was derived using 
stars between 9° — 13°, provides quite a good fit to the den¬ 
sity profile over the full radial range out to 18.5 kpc (~ 20°) 
despite the apparent presence of the faint western structure. 

To investigate the region to the west of £ = 0 in more 
detail, we constructed field subtracted Hess diagrams in the 
radial range 12.5 — 14.5 kpc and 14.5 — 18.5 kpc. We utilised 
exactly the same procedure as that which we previously used 
to construct the field subtracted Hess diagrams for the main 
body of the stream-like overdensity (Figure [5} , but in this 
case adopting field boxes of appropriate area from the west¬ 
ern region just outside 18.5 kpc. These Hess diagrams are 
presented in Figure [13] and reveal that LMC-like popula¬ 
tions are clearly visible even to very large radii as inferred 
from the stellar density map and radial profile. To the preci¬ 
sion possible with such low-density CMDs, the populations 
in this region appear to match both the LMC disk within 
12.5 kpc (« 13.5°) and the main stream-like overdensity. 

Returning to the stellar density map and now consider¬ 
ing the region to the east of £ = 0, the picture is complicated 
by the presence of foreground contamination that increases 
noticeably towards the edge of the survey footprint. This is 
due to the Milky Way disk, which lies roughly parallel to 
the edge of the survey footprint another ~ 20° further east. 
We have already discussed the issue of foreground contami¬ 
nation in Section 13. 3 1 here, we note that these contaminants 
make a non-negligible contribution to the surface density 
profile outside ~ 16.5 kpc (marked with a grey dashed el¬ 
lipse on the map and a vertical dotted line on the profile), 
which is the radius within which the main stream-like fea¬ 
tures are concentrated east of £ = 0. The blue profile in 
Figure ITTI should not be trusted outside 16.5 kpc; this con- 
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Figure 13. Field subtracted Hess diagrams for stars lying to the 
west off = 0 and between 12.5—14.5 kpc (left) and 14.5—18.5 kpc 
(right). The fiducial isochrone and CMD selection boxes are the 
same as in Figures [I] and [5] LMC stellar populations are evident 
even to very large radii. 


tamination is the reason why the profile for all stars, in the 
upper panel of the plot, does not converge to the background 
level (determined using the profile to the west of £ = 0) until 
very large radii. 

In summary, we find that to the west of £ = 0 the 
LMC radial surface density profile shows clear evidence for 
very diffuse populations reaching out to 18.5 kpc (« 20°). 
These may also be present to the east of £ = 0, but we 
cannot be certain due to increasing contamination from the 
Galactic foreground. Our observation of such remote LMC 
populations is consistent with the spectroscopic detection 
in many dire ctions of prob a ble LMC g iants out to radii 
of ~ 2 0° by iMaiewski et al.1 (|200(A l2009h and iMunoz et al.1 
( 200( 1), and the recent blue- horizontal branch star maps of 
iBelokurov fe KoposovI ( 2016 ). Our overall LMC surface den¬ 
sity profile shows a gentle break outside a radius of 12.5 kpc 
from the best-fit exponential fall-off measured in Section [331 
to a flatter decline. T his is quite similar to the profile ob¬ 
served by IMaiewski et ahl (|2009) - although they assume an 
exponential scale length of 1.6° ~ 1.5 kpc and a break radius 
of 9° ~ 8.5 kpc, their spectroscopic data points (top panel of 
their Figure 2) exhibit a substantial spread and would also 
be compatible with our slightly longer scale length of 1.64 
kpc and larger break radius of 12.5 kpc ~ 13.5°. In addition 
we find that these remote populations have characteristics 
similar to those of stars in the outer LMC disk and in the 
stream-like overdensity we have discovered stretching to the 
east. This observation i s again compatible wit h the spectro¬ 
scopic observations bv IMaiewski et all ( 2009 ). which show 
an apparently flat metallicity profile out to ~ 20°, with a 
mean value of [Fe/H] ~ —1 .0. _ 

While Maie wski et al.l ([20091 ) interpret the break in 
their profile to indicate the presence of a classical halo 
around the LMC, the break in our overall profile can be 
fully explained by the stream-like feature extending east - 


when we construct a profile using only stars to the west 
of f = 0, where this overdensity is not present, it is well 
described by our best-fit exponential disk model over the 
full radial range out to 18.5 kpc (« 20°). This does not 
rule out the presence of a stellar halo - it is important 
to note that the fields observed by IMaiewski et al.l ( 20091 ) 
span a signfiicantly wider range in azimuth than does the 
present D ES footprint , and, moreover, the recent maps by 
IBelokurov &; KoposovI ( 2016 ) tracing blue horizontal branch 
stars suggest an LMC halo component almost certainly ex¬ 
ists. However, it does indicate that (i) the LMC exponential 
disk may extend further than previously suspected, at least 
in some directions, and/or (ii) the transition to a halo-like 
population is rather smooth such that it does not lead to a 
noticeable break in the profile (at least to the north-west). 
If the disk indeed extends out to such large radii, this would 
be in contrast to some previous estimates for its tru ncation 
radius, which he in the ran ge ~ 13 — 15 kpc (e.g. JSaha et al.l 
l20ld : iBalbinot et al.l 120151 ) . We have already noted in Sec¬ 
tion GU how the sharp drop off in density ~ 15° — 16° due 
north of the LMC cou ld e xplain the apparent trunction in 
the profile of I Saha et al.l (2010). and we suspect a similar 
effect could explai n the “edge” seen at ~ 13.5 ± 0.8 kpc to 
the north-west by IBalbinot et al.l (|2015l ). In this direction 
we observe the faint western structure extending to ~ 14.5 
kpc; it may be that our CMD selection box, which traces 
much more of the main sequ ence a nd main sequence turn¬ 
off region than that of IBalbinot et all (120151) . has facilitated 
the detection of extremely diffuse populations outside this 
radius - note that on the CMD in Figure fl3l LMC stars be¬ 
yond 14.5 kpc are only clearly detected below the turn-off. 


4 DISCUSSION & SUMMARY 


We have used year one DES imaging to study the structure 
of the northern periphery of the LMC. We find evidence 
for a substantial stream-like stellar overdensity stretching 
from ~ 13.5° due north of the LMC centre towards the 
east. Although this is, broadly speaking, a similar direc¬ 
tion to that of the leading arm of the Magellanic Hi stream, 
there is no obvious correspondence between the two - the 
gase ous feature is much w ider and more extended (see e.g., 
iNidever et al.l 120081 . l20ld ) , and there is at present no sug¬ 
gestion that our substructure comprises a stellar component 
to it. 

The width of our overdensity is ~ 2°, or w 1.5 kpc, 
and it spans an arc on the sky that is, at minimum, 12.5° 
long - corresponding to a physical length greater than ~ 10 
kpc. This assumes that the feature passes in front of the 
Carina dwarf galaxy around 20° north-east of the LMC, 
possibly explaining why a number of previous studies of this 
object have found evidence for LMC po pulations in their 
field of view (e.g .. iMcMonigal et al.ll20li ) that exhibit cold 
kinematics (e.g., IMunoz et al.l 120061 ). We can think of three 
possible scenarios for the origin of this structure: 


• It might be the remains of a disrupted satellite of the 
LMC. Recent searches have uncovered a substantial number 
of low-luminosity dwarf galaxies in the vicinity of the LMC 
jKojDOsa^e^jd. 2015; Bec htol et al. 20151 : iKim et al.l 120151 : 
Kim & Jerien 2015; Martin et al. 2015), some of which may 
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well be bound to the Magellanic Clouds le.g.. iDeason et ahl 
120151) . 

• It may be an overdensity in the remote outskirts of the 
LMC disk itself - perhaps akin to an extended s piral arm. 
The recent diffuse light imaging by iBesla et~ahl (2016) has 
revealed similar structures extending to « 8° from the LMC 
centre. 

• It could be comprised of material that is being stripped 
from the outer LMC disk due to the tidal force of the Milky 
Way, or due to a recent close encounter with the SMC. 


The observed properties of the substructure help shed 
more light on these possibilities. We find that the integrated 
luminosity of its brightest portion is My ~ —7.4. This sits 
in between the Hercules dwarf (My ~ —6.6) and the Draco 
and Ursa Min or dwarfs ( My ~ —8.8) in the Milky Way sys¬ 
tem (e.g., iMcConnachiel [2012 1. although we note that our 
measurement is a lower limit since we have not traced the 
full extent of the overdensity on the sky. Our field-subtracted 
Hess diagrams show that the stellar populations in the sub¬ 
structure are predominantly old, as would be expected for 
a faint dwarf, but, critically, they also have [Fe/H] ~ —1. 
This is much more metal-rich than would be expected for 
such a low luminosity syste m acco rding to the luminosity- 
metallicity relation (e.g.. iKirbv et al1l201lh - Draco, Ursa 
Minor, and Hercules all have [Fe/H] < —2. It could be that 
the putative satellite was initially much more massive and 
has suffered extensive tidal destruction; however the original 
system would need to have been unfeasibly large - nearly 
two orders of magnitude more luminous than at present, 
comparable to the Fornax dwarf with My ~ —13.4. 

Our photometry reveals that, in fact, the stellar popu¬ 
lations seen in the stream-like feature are an excellent match 
for those found in the outer LMC disk. This is strongly sug¬ 
gestive that it is comprised of disk material - perhaps stars 
being tidally stripped from the disk, or an overdensity in the 
structure of the disk itself. Many previous studies have found 
non-axisymmetric features in density maps of the LMC - 
there are spur-like struct ures, possibly spiral arms, in the 
inner regions (e.g.. IZaritskv et al.1 12004 1: while at interme¬ 
diate radii t here are spur-like and shell-like features present 
(e-g., i rwi ami), some of whic h sit t o wards the north and 
progr ess from west to east (e.g., de Vaucouleurs &; Freemanl 
119721 1. More recently, IBesla et all (|2016h have revealed arc- 
like structures in the disk extending to ~ 8°. Furthermore, 
we clearly observe additional structure of some description 
in the LMC disk at radii of 9° —12° (see Figure flQl). although 
the nature of these features is not clear at present. Despite 
the fact that we detect no young populations at the radius of 
our overdensity, such that it is probably not a spiral arm in 
the traditional sense, it is not difficult to imagine a scenario 
in which one or m ore close en counters between the LMC and 
SMC (as in IBesla et al.ll20l'6 U or between the LMC and the 
Milky Way, could give rise to transient “ripples” or pertur¬ 
bations in the outer LMC disk, or indeed remove material 
entirely into a stellar stream. 

In either of these cases, it might be expected that the 
outer LMC disk ought to show some other signs of distor¬ 
tion. We have observed tentative evidence for this. By con¬ 
structing models of inclined planar disks we find that the 
observed stellar density of the outer LMC disk on the sky, 
and its azimuthal distance profile, can only be simultane¬ 


ously explained by an elongated disk. The necessary ellip- 
ticity is only mild, but might in dicate the presence o f tidal 
distortion (see the discussion in Ivan der Marelll200lh . It is 
also relevant that the orientation of the elliptical isodensity 
contours of the outer LMC is such that the major axis is 
aligned very close to the north-south direction on the sky, 
coinciding with the apparent location of the origin of the 
stream due north of the LMC centre. 

In addition to this, the position angle of the line of nodes 
in our best-fitting elliptical disk model is twisted relative 
to all previous determinations from more central regions, 
which may indicate a warp in the disk. Our measured in¬ 
clination is also lower than most estimates derived from the 
inner LMC; indeed, it follows the general trend of decreasing 
i nclina tion with radius observed by Ivan der Marel fe Cionil 
(l200lh . Finally, there is an offset in our absolute distance 
measurement to the outer LMC disk, which we find to be 
about 0.15 — 0.20 mag brighter than predicted by our best¬ 
fitting planar model. While much of this offset can likely be 
attributed to systematics associated with deriving distances 
using isochrone models, it might also reflect warping or flar¬ 
ing of the disk, or some other kind of distortion along the line 
of sight. This type of non-planar behaviour at large radial 
distances is again possibly indicative of tidal perturbations 
or distortions in the outer disk. 

If our substructure is indeed comprised of stripped disk 
material, then its origin might be due to the tidal force 
of t he M ilky Way. This is the argument put forward by 
Ivan der Marell (|200ll l to explain the (substantial) intrinsic 
ellipticity they infer for the LMC disk - the orientation that 
they measure for the elongation points within ~ 6° of the 
direction of the Galactic centre, which sits at a position an¬ 
gle of 183.7° relative to the LMC centre. Our results are in 
excellent agreement with this finding - the position angle of 
the line of nodes in our best elliptical disk model sits at 185°, 
while the ellipticity on the sky is almost perfectly aligned in 
the north-south direction (see e.g., the lower panel of Figure 
E}. Moreover, as expected in the case where our substruc¬ 
ture is a stream induced by the tidal force of the Milky Way, 
its direction is roughly perpendicular to the disk elongation 
and quite well aligned with the latest Hubble Space Telescope 
proper motions, which are fiw — —1.910 d= 0.020 mas yr" 1 
and fiN = 0.229 d= 0.047masyr -1 ([Kallivavalil et all l2013h 
- i.e., implying a tangential motion predominantly towards 
the east. Observations from the Pan-Andromeda Archaeo¬ 
logical Survey (PAndAS) have revealed the presence of sim- 
i lar stellar streams emana ting from the di sks of both M33 
(iMcConnachie et ~aT1 l2010h and NGC 147 ([Crnoievic et al.l 
|2Qi4 see also Irwin et al. 2016, in prep.). These are inferred 
to have arisen due to recen t encounters between thes e satel¬ 
lite galaxies and M31 (e.g. iMcConnachie et al.ll2009h . 

To explore this possibility in more detail we ran sev¬ 
eral iV-body simulations of an LMC flyby. These were car¬ 
ried out using the iV-bod y com ponent of gadget -3 which 
is similar to GADGET -2 (lSpringelH2005t ). In each simulation, 
the LMC was modelled as a two com p onen t galaxy w ith a 
stellar disk and an NFW profile ([Navarro et al.|[l997l b We 
constructed self-consistent and s table initial conditions fol - 
lowing the procedure described in Villalo bos k, Helmil {2008) 
and used the same form for the stellar disk as in that work. 
Before adiabatic contraction, the LMC halo had a mass of 
M200 -1.2 x 10 11 Mq and a concentration of C200 — 8 . The 
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stellar disk had a mass of Md = 4 x 10 9 Mq, a scale ra¬ 
dius of 1.5 kpc (matching the inferred value in Table 1), 
and a scale height of 0.3 kpc. After the adiabatic contrac¬ 
tion of the LMC halo, the total LMC mass within 8.7 kpc 
is 1.8 x 10 10 Mo, consi stent with the mass inferre d fr om th e 
LMC rotation curve in Ivan der Marel fe Kallivavalill (120141 . 
The total LMC mass is 1.4 x lO 11 !!^, consistent with 
other recent mass esti mates (e.g.. IKallivavalil et al.l [20131 : 
IPenarrubia et all 120151 ) . The initial orientation and rota- 
tional sense of the LMC were chose n to match the results of 
Ivan der Marel &; Kallivavahl (|2014h . We used 10 6 particles 
each to represent the LMC disk and halo, and a softening 
length of 75 pc and 500 pc respectively. The Milky Way was 
modelled as a 3-component system with a_bulge, disk, and 
dark matter halo as described in I Gomez et al.l ( 20151 ). 

Our iV-body model was integrated for 2 Gyr. The initial 
positions and velocities of the Milky Way and LMC were 
chos en using backward integration from the current position 
as in I Gomez et al 1 (I 2015 I') . which resulted in the LMC being 
wi thin 2 a of the present d ay positions and velocities given 
in IKallivavalil etahl ( 20131 ). 

Figure [2] shows the density of particles in the disk - 
i.e., the stellar density in the simulation - using the same 
coordinate system as in Figure [2] The simulation shows ma¬ 
terial being stripped from the northern part of the disk in 
the correct direction to give LMC debris near the location of 
the Carina dwarf - thus at least qualitatively reproducing 
the properties of our stellar overdensity. In the model the 
stream debris is due to tidal stripping of the LMC, and the 
rate of this stripping is controlled by both the orientation 
of the LMC inte r nal rotation and its o rbital motion (e.g., 

I Read et all 120061 : iD’Onghia et al.l l2009h . As a consistency 
check we re-ran the simulation with the LMC in two orien¬ 
tations, prograde and retograde, and found enhanced and 
decreased stripping respectively, as expected. For the orien¬ 
tation used in our simulation, the direction of the LMC spin 
and the angular momentum of its orbital motion around the 
Milky Way differ by 75°, so we expect a slight enhancement 
of the stripping over a non-rotating progenitor. Intriguingly, 
at the time shown in Figure fTH our simulated LMC disk is 
intrinsically elliptical, with an ellipticity e ~ 0.15 at radii 
commensurate with those considered by our disk-fitting al¬ 
gorithm in Section 1331 

We emphasise that these simulations do not prove that 
the observed substructure is a debris stream caused by tidal 
interaction between the LMC and the Milky Way, merely 
that this is a plausible explanation. In the future we plan 
to run a larger suite of models to explore how the debris is 
affected by varying the properties of the Milky Way and the 
LMC, and the presence of the SMC. 

The SMC is certainly a complicating factor, and indeed 
we cannot rule out that our substructure is due to one or 
more close encounters between the two Magellanic Clouds, 
rather than simply due to the tidal force of the Milky Way. 
The SMC sits at a position angle of 228.9° relative to the 
LMC centre, which is broadly in alignment with the direc¬ 
tion of the arc of our stream-like overdensity. This galaxy 
is believed to have had multiple close interactions with the 
LMC in the past - perhaps as recently as a few hundred Myr 
ago, and perhaps passing as close as a few kpc from the 
LMC centre. Observations have revealed a stellar co mpo- 
nent in the Magellanic Bridge region fe.g.. llrwin et al1ll985l : 
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Figure 14. Stellar density of our simulated LMC using the same 
coordinate system as in Figure [2] The northern outskirts of the 
disk show a stream of tidally stripped debris qualitatively similar 
to our observed substructure. 

Nidever et al.l 120131 : iNoel et all 120131 . 120151 : ISkowron et al.l 
20141 ). probably stripped from the SMC by just such an en¬ 
counter. 

iBesla et"ahl (l2012l j model the effects of a very close 
LMC-SMC interaction, where the orbits are consistent with 
the observed proper motions; amongst their findings is that 
such an interaction can substantially disturb the LMC disk 
and even give rise to possible stream-like features - in the 
particular case they consider this mainly occurs to the north 
and north-eas t of t he LMC cen tre as we have observed. 
More recently, IBesla et"al1 ( 20161 ) have conducted modelling 
in light of their discovery of stellar arcs and spiral arms 
extending to ~ 8° north of the LMC, and the absence of 
comparable features in the south. They show that repeated 
close passages of the SMC can lead to precisely this type 
of asymmetric arrangement. Together these computations 
show that, at least in principle, our LMC overdensity could 
arise due to interaction(s) with the SMC; moreover, such en- 
couters might also offer an explanation for the structure we 
observe in the LMC disk between 9° — 12°. Ultimately, new 
observations and more detailed modelling will be required 
to identify the most likely origin for our stream-like feature. 
A key question is whether a matching substructure exists 
due south of the LMC centre and stretching to the west - 
the tidal features seen around both M33 and NGC 147 due 
to the influence of M31 are q ui te sym metric in this regard, 
while the models of lBesla et al.l ([2016 9 show that if the SMC 
is responsible then the arrangement should be asymmetric 
(i.e., there should be no matching arc in the south). 

We note that radial velocity observations along our sub¬ 
structure should help determine whether this feature is a 
stream of material stripped from the outer LMC disk, or an 
overdensity in the disk itself. If the latter, then it ought to 
approximately follow the known disk rotation profile with 
position angle. If, on the other hand, it is a stream of tidal 
debris, then our iV-body model predicts a much more con¬ 
stant velocity with position because the stripped material 
tends to simply proceed roughly along the LMC orbit rather 
than continuing to share in the disk rotation pattern. 
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Finally, we have observed apparent evidence for LMC 
populations at radii of up to ~ 18.5 kpc (« 20°) to the 
north and north-west. This corroborates the results of pre¬ 
vious studies that have spectroscopically detected apparent 
LMC populations in discrete fields at radii of u p to ~ 20° 
and over a wide range of azim uthal angles fe.g.. lMuhoz et all 
2006; iMaiewski et, all 12009 'll we are now beginning to cre¬ 
ate a fully-filled map of these regions. The nat ure of these 
extrem ely remote stars is not clear - while IMaiewski et al.l 
(2009) interpret them as a classical spheroidal halo, our ra¬ 
dial density profile to the north-west is well described by our 
best-fitting inclined disk model over all radii beyond ~ 9° 
without evidently requiring a separate halo component. Fur¬ 
thermore, it is not obvious whether these extremely remote 
stars could still be gravitationally bound to the LMC - we 
have observed material that may very well be in the process 
of being tidally stripped from the LMC at a radius of ~ 12.5 
kpc from the galaxy’s centr e; _on the other hand , the kine¬ 
matic measurements of Ivan der Marel fc Kallivavalill (l2014h 
suggest that the tidal r adius o f t he LMC cou ld be as large 
as 22.3 d= 5.2 kpc, and iBelokurov &; KoposovI (2016) have 
observed possible LMC members extending to even larger 
radii. Whether or not these outermost populations remain 
bound to the LMC is a particularly interesting question - 
if they are not, then the fact that we are still able to ob¬ 
serve their presence would likely argue for the LMC only 
having entered the tidal influence of the Milky Way quite 
recently, such that it is possibly still on its first pass around 
the Galaxy. Further deep imaging and spectrscopy is neces¬ 
sary to explore the properties of this very diffuse structure. 
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